How does transverse (hydrodynamic) flow affect jet-broadening and jet-quenching ? 



R. Baier 1 *, A. H. Mueller 2 t and D. Schiff 3 * 
1 Physics Department, University of Bielefeld, D-33501 Bielefeld, Germany 
2 Physics Department, Columbia University, New York, NY 10027, USA 
3 LPT, Bat. 210, Unwersite Paris-Sud, F-91405 Orsay, France 

February 9, 2008 



^ ■ We give the modification of formulas for pj_-broadening and energy loss which are necessary 

' to calculate parton interactions in a medium with flow. Arguments are presented leading to the 

£nQ , conclusion that for large px-spectra observed in heavy ion collisions at RHIC, the influence of 

transverse flow on the determination of the "quenching power" of the produced medium is small. 
This leaves open the question of the interpretation of data in a consistent perturbative framework. 

S: 

, 1. Introduction 

■ 

Energy loss of a high transverse momentum parton travelling in the hot and dense medium created in ultrarel- 
ativistic collisions has been the subject in recent years of intense investigation, reviewed in [1-6]. With the aim of 
\Q , being able in particular to assert the existence of the QGP phase of matter [7] . 

It has become clear however, that an important issue is the influence of the medium evolution on the radiative 
energy loss, when following the BDMPS [8-10] - Zakharov [11,12] - Wiedemann [13] approach. The medium cannot 
be described as being static [14] : longitudinal and transverse flow have to be considered [15-17]. They may alter 
substantially the distribution of matter before freeze out. The crucial parameter for energy loss is the local transport 
coefficient q{r) , related to the squared average transverse momentum transfer - from the medium to the hard parton 
~h | - per unit length. 

For an ideal QGP, and in a comoving coordinate system, one may relate locally q with the energy density of the 
Q ' medium e : q ~ c e 3 / 4 with c ~ 2 [18, 19]. 

Recently from various authors [20-23] came the determination of c from data, taking into account Bjorken lon- 
Ch ■ gitudinal expansion [24] but neglecting transverse flow. In [20] c is found to be much larger than 2 (c > 8 • ■ • 19, 
c ~ 10), leading to speculations about a strongly coupled QGP [25-30] . A recent endeavour [31,32] implementing a 
strong non- Bjorken expansion including radial flow [33], with a small value of primordial transverse velocity v % = 0.1 
has led to a "moderately optimistic" scenario reducing c = 10 to c = 2 compatible with perturbative estimates. 
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In the present note, we try to reformulate the problem. We show first how to determine q in the presence of 
transverse flow by applying a proper Lorentz boost. We assume a realistic (ideal) hydrodynamical description for 
heavy-ion collisions, based on longitudinal Bjorken expansion [24] and a radial flow with vanishing initial velocity. 

We are led to the conclusion that the transverse flow is a small effect and does not provide the solution to the 
perturbative /non-perturbative dilemma. 



2. Moving medium versus medium at rest : geometry 



The problem is to find the relationship between <f for a moving medium and go (medium at rest), locally in space 
and time. We focus on central collisions. 
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We consider - in the transverse plane (xy) of a collision- a medium moving with velocity v along the y axis towards 
a parton which enters into it at the origin of our coordinate system at time r = 0, at point A, leaving it at point B 
with = i(l,sin0, cos0, 0) where we denote by t the time spent by the parton in the medium. The dimension of 
the medium in the velocity direction being L, one finds 

t = ^— • (2.1) 

v + cos 

The medium at rest has length = = L ch £ defining £ such that v = th (. Entering at point A , the parton 

VI — v 2 

leaves the medium at rest, at point B such that as a result of the Lorentz boost, we find 

Bq = t(ch C(l + cos0 th C),sin0, ch C(cos0 + th £),0) (2.2) 
with 

AqBq = L = t ch C (1 + cos0 th C) (2.3) 
cos 9 a 

The angle 9 corresponds to the angle 9, but evaluated in the rest system of the medium, implying 

. th ( + cos0 . . 

cos 6» = — „ : - (2.4) 

1 + cos 9 th C 



and 



Sm9 ° = chC(l + cos^thC) • (2 ' 5) 



3. Transverse momentum broadening 
A. Local geometry 

Let us define in the rest frame the initial momentum of the parton : 

Po = Po(l, sin O , cos 9 ,0) . (3.1) 

After travelling on a small distance St - we want to extract a local information - and experiencing p± -broadening 
and energy loss, the parton has a momentum : 

Pq = (Po - Seo) (l,sin (0 O + <50 o ) cos^o, cos (0 O + 59q) , sin (0 O + S9 ) sin5<p ) (3.2) 
where 5ipo is the azimuthal angle in the zx plane. In the moving frame, the corresponding momenta are : 

p» =p(l,sin0,cos0,O) (3.3) 
p» = - fe) (1, sin (0 + 50) cos ^, cos (0 + ^0), sin (0 + ^0) sin ^) . (3.4) 

The momenta in both frames are related via a Lorentz boost, which leaves the x component unchanged, leading to 

psin0 = posin0o • (3-5) 

On the other hand dip = S(po, since the Lorentz boost does not affect the projection on the zx plane. Using eqs. 
(2.4) and (2.5), we also find : 

69 590 (3.6) 



sin 9 sin 9 a 
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B. Comparing frames 



The corresponding transverse momentum broadening is easily obtained in both frames : we sum the squares of 
the momentum increase in the z axis direction and in the direction orthogonal to the parton momentum : 

5pl± = Sp 2 Qz + Sp' 2 Q = pi [sin 2 Q Sip 2 + 59 2 Q ] . (3.7) 
Similarly 

Spl = p 2 [sin 2 9 Sip 2 + 69 2 ] . (3.8) 

Using eqs. (3.5) and (3.6), we find that Sp\ = 5p^ ± and may immediately infer the relation between q and q : indeed 
we get 

5p\ =q5r and (3.9) 
Spl± = 90 St = q St [ch C (1 + cos 9 th ()] , (3.10) 

leading to (locally) : 

q = q ch C(l + cos^thC) . (3.11) 



Recently this transformation property has been independently derived in [34], and already used in [35]. 
§b may be parametrized as discussed in the Introduction. 

Writing the transverse flow velocity as u M = (ch (, 0, — sh (, 0) and with the parton momentum written as 

Pa 

( 1, sin 9, cos 6,0), we get 

5 = 50 — . (3.12) 

Po 

A way to understand this relation is to write the transport coefficient as the rate of interaction - the number of 
interactions per unit time - multiplied by the typical invariant scale /i 2 characterizing the momentum transfer to the 
parton in a parton-medium collision : for a medium at rest % = Rq[i 2 and for a moving medium q = R/j, 2 . The rate 
of interaction varies indeed like the inverse of the time spent in the medium : 

■J- = ch C(l + cos 6> th C) , (3.13) 
Ro 

so that eq. (3.11) is recovered. 



C. Parton transverse broadening 



Using qo(r) we may now calculate the transverse broadening [9] obtained when taking into account the transverse 
flow for a parton travelling in a thcrmalizcd medium at temperature T. This last condition allowing us to determine 
the r dependence of T and the transverse velocity v by using hydrodynamical model codes. 

In Fig. 1 we have schematized the situation in the transverse plane [20] : the vector s(s, <fi) defines the point in the 
plane where the hard parton enters the medium. The parton then follows a path along the unit vector n and along 
its way, encounters the radial flow along the direction r which depends on time r. 
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Fig. 1. View of the kinematics in the transverse plane. 

We write f = s + rn and denning 9 (according to the definition given in sect. 2) and <f> as in Fig. 1, we find that 

s cos 6 — t 



cos 9 



Writing ch £ as 



S 2 + T 2 — 2STC0S 
1 



(3.14) 



j(v) and th £ = v, where v depends on r and r. We may now get the following expression 



for the transverse broadening of the parton for a cylindrical medium at central rapidity : 

(Ap^Bj+fiow = sr I sds d< t> I dT *) (T(r,T)) j(v) (1 + v cos 9) 
nn A J J J To 



(3.15) 



where T max is defined as the maximum in-medium path length compatible with geometry and taking into account 
the life-time of the thermalized medium tqgp ■ 



r m ax = min{r^ e a " i ,TQ GP } 



with 



s 2 cos 2 6 + R 2 , 



(3.16) 



(3.17) 



and To is the initial time. In the integral (3.15) the boundaries are simplified by neglecting the flow beyond the 
nucleus radius Ra (c.f. Fig. 3). The life-time tqgp is estimated in the framework of Bjorken expansion [24], after 
having fixed the initial temperature and taking the final temperature as T ~ 200 MeV. (In the following we take 
tqgp — 4.5 fm). 

We actually calculate the ratio 



Rflow = (ApjjBj+flow/i^P^Bj 



(3.18) 



in order to quantify the importance of the radial flow in addition to Bjorken (Bj) longitudinal expansion [24], for 
which 



4o(T) = «o(T )(-) = ( fo(r )(^) 3c = 

JO T 



(3.19) 



where c s denotes the sound velocity (taken to be c s = -^=). 

The calculation is done in the framework of ideal hydrodynamics, assuming initial conditions relevant for heavy-ion 
collisions, e.g. a vanishing flow velocity at initial time t [36]. The result is shown as the solid curve in Fig. 2 for 
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a cylindrical medium of radius Ra = 6.4 fm {Au nucleus), fixing the initial time tq = 0.5 fm, as a function of the 
initial temperature T . 

The curves are actually calculated using the approximate analytic expressions derived in [36] (see also [37]). In 
the ratio Rfi ow the dependence on the initial temperature To cancels. Ignoring the limit given by tqgp results in 
the dotted curve. We note that Rfi ow is smaller than 1. This can be understood in the following way: the flow has 
a non negligible effect for large enough values of r, where v differs significantly from 0. But this is only realized (see 
Fig. 1) when the jet is moving with the flow, since then cos 9 ~ —1, and therefore j(v)(l + vcos9) < 1. This effect 
is even larger for the dotted curve. 

We note that the angular dependence of eq. (3.11), which is responsible for the reduction of broadening, and also 
quenching, is not present in the ansatz for q given in [33] and used in [31,32]. 



l.i- 



l - 



9 



0.8 - 

, i , 

4 6 8 

T [GeV] 

Fig. 2. Rfiow as a function of initial temperature To from eqs. (3.15) and (3.18) (solid curve). Dashed and dotted curves 
as described in the text. 

As an exercise we calculated Rfi ow with the following set-up: we assume that qh(T, r) follows Bjorken's r de- 
pendence (3.19), keeping the hydro velocity v(t,t) in 7(w), but without the r-dependence in the temperature, and 
dropping the cos 9 term. The ratio Rfi ow becomes larger than one, which can be seen from the dashed curve in 
Fig. 2, because in this case the factor j(v) wins. 

Two conclusions emerge : 

a) the effect of the transverse flow is small and of the order < 10 %. 

b) the ratio Rfi ow is smaller than 1 which goes against the possibility that the transverse flow would allow us to 
solve the perturbative/non-perturbative dilemma discussed in the introduction. 

The smallness of the effect is due to the fact that the transverse velocity takes a long time to reach an appreciable 
value, starting from v = and cannot provide the possibility of compensating the cooling effect. This can be 
visualized in Fig. 3 where we have plotted the quantity (T/T)) 3 • j(v) - which is essentially the integrand in eq. 
(3.15) - for different values of r, as a function of r, using ideal hydrodynamics [36]. 

It is interesting to investigate the effects due to a non-vanishing shear viscosity compared to the ideal hydrody- 
namical description for the quantity (T/T)) 3 • 7(w), which is plotted in Fig. 4. The details of this calculation may be 
found in [37]. 
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Fig. 3. (T/Tq) 3 • 7(«) as a function of r for different values of r calculated from ideal hydrodynamics [36] 
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Fig. 4. Effects due to non- vanishing shear viscosity rj [38]. 

No significant effects due to the transverse flow in the presence of shear viscosity are found, as maybe expected, 
when comparing the solid and dashed curves in Fig. 4, although cooling is slowed down at least for times r presented 
in this figure: Rfi ow stays of 0(1). 
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4. Energy loss 



The medium-induced radiative jet energy loss is determined by the transport coefficient q. While the energy loss 
formalism has been previously developed only for a medium without flow it is apparent from the discussion at the 
end of section 3.B that if one uses q, as given in eq. (3.11), to determine the interaction of QCD partons with the 
medium then there is no modification of the fundamental energy loss formulas. Thus eqs. (25) and (26) of [15] remain 
valid in the presence of flow. The basic quantity is the radiation spectrum of gluons emitted from the high energy 
parton. In the soft gluon energy limit, to much smaller than the parton energy, the dominant mechanism consists in 
multiple scatterings of the radiated gluon off the "centers" describing the medium [8]. For static centers, as well as 
for a dense (thermal) medium which undergoes longitudinal Bjorken-type expansion [24] the gluon spectrum obeys 
the following scaling law in terms of a characteristic gluon energy lo c [2, 16], 

W £ = I(u/w c ) . (4.1) 
For static centers [8] 

= \qL 2 , (4.2) 

where L denotes the path length of the energetic jet in the medium. 
In the expanding case [15], eq. (4.2) is generalized to [16,17] 



Wc= / dTTq{r) , (4.3) 
Jo 

with q(r) s=s 1/r, when following Bjorken [24]. 

However, when transverse flow is present in the medium this nice scaling property does not hold, since q appears 
at different times, corresponding to the interference product of the emission amplitude and the complex conjugate 
one (see eg. eq. (25) in [15]). As for the case of p± -broadening, we do not expect large effects due to transverse flow. 
Therefore, for simplicity, we calculate the ratio 

^ = (4 4) 

In comparison with the calculation of Rfi ow , we only estimate the first moment with respect to the path length r of 
the energetic jet, i.e. 

(u c )Bj+fiow = sr sds d( t> I dT T qo(T(r,T))"f(v) (1 + v cos 6) , (4.5) 
Ttlt A Jo Jo Jo 

and accordingly for (u> c )Bj- 

As we expect from Fig. 2 the ratio R UJc does not differ significantly from 1, when calculated within (ideal) 
hydrodynamics and realistic initial conditions for the flow field v(r,To): a typical value is Ru, c ~ 0.85 for 
t = 0.5 fm, To = 400 MeV. From this estimate, we expect small effects of radial flow on the quenching of 
large transverse momentum hadrons produced in nucleus-nucleus collisions. This conclusion on R Uc < 1 is confirmed 
in the recent analysis in [35] using eq. (3.11). 

Our study seems to rule out the possibility of reducing the large value of q from model comparisons with data [20] 
by including radial flow with initial conditions, which we assume to be realistic and causal in heavy-ion collisions. 
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